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Near-field interactions and non-universality in speckle patterns 
produced by a point source in a disordered medium 
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A point source in a disordered scattering medium generates a speckle pattern with non-universal 
features, giving rise to the so-called Co correlation. We analyze theoretically the relationship between 
the Co correlation and the statistical fluctuations of the local density of states, based on simple 
arguments of energy conservation. This derivation leads to a clear physical interpretation of the Co 
correlation. Using exact numerical simulations, we show that Co is essentially a correlation resulting 
from near-field interactions. These interactions are responsible for the non- universality of Co, that 
confers to this correlation a huge potential for sensing and imaging at the subwavelength scale in 
complex media. 
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Waves propagating in a disordered medium generate 
a strongly fluctuating spatial distribution of intensity 
known as a speckle pattern [l|. The spatial structure 
of a speckle pattern is often characterized by the inten- 
sity spatial correlation function (/(r)/(r')), or its an- 
gular counterpart (/(u)/(u')) where the unit vector u 
defines an observation direction. In usual experiments 
the medium is illuminated by an external beam, and the 
speckle pattern is observed, e.g., in transmission. Short- 
range and long-range contributions can be identified in 
the intensity correlation function, which is written as a 
sum of three terms denoted by Ci (short range), C2 and 
C3 (long range) [2] . These correlations have been widely 
studied since they are responsible for enhanced meso- 
scopic fluctuations [3| and their sensitivity to changes 
in the medium can be used for imaging in complex me- 
dia [4|. Moreover, the possibility of controlling speckle 
patterns by wavefront shaping has generated new inter- 
est [5| . When the waves are generated by a point source 
placed inside the medium, a new type of spatial corre- 
lation appear, that has been denoted by Co iBj. This 
correlation has features that make it particularly inter- 
esting: It is of infinite range and can dominate the usual 
long-range correlations, and it is strongly dependent on 
the local environment of the source [y, |7|. This non- 
universality makes Cq a valuable quantity for sensing or 
imaging in complex media with a high sensitivity to the 
microscopic structure of the medium [7|. Moreover, it has 
been shown that the Cq correlation and the fluctuations 
of the local density of states (LDOS) at the location of 
the point source are equal [8| . This means that Cq could 
be obtained from measurements of LDOS fluctuations in- 
stead of an analysis of speckle patterns. In optics, LDOS 
fluctuations can be measured from the fluorescence life- 
time of nanoscopic emitters |9l4l3l|. They are expected 
to provide information on the local environment of the 
emitter [ij, [l5| and on the photon transport regime |16l - 
|20|. In this Letter, we revisit the Cq correlation concept 
and its connection to LDOS fluctuations in the context 



of scattering of electromagnetic waves. We consider elec- 
tromagnetic waves since the most recent relevant experi- 
ments have been performed in this regime. First, we give 
a novel derivation of the equality between Co and the nor- 
malized variance of the LDOS based on energy conserva- 
tion, showing that the relationship is valid in any regime, 
including strong localization 1] . This derivation leads for 
the first time to a clear physical interpretation of the Co 
correlation. Second, we analyze the non-universality of 
Co, and show that it is due to a large extent to near- 
field interactions. The sensitivity to the degree of cor- 
relation of disorder around the location of the source is 
demonstrated based on numerical simulations on three- 
dimensional open systems. 

In this first part, we define the Co correlation and 
derive its relationship with LDOS fiuctuations based 
on energy conservation. We consider a non-absorbing 
disordered system of finite size, embedded in a sphere 
with radius R (this geometry is that used in the subse- 
quent numerical simulations). The system is illuminat- 
ing by a classical point dipole p radiating at frequency 
cj, placed at the center of the cluster at position r^. The 
medium is assumed statistically isotropic and homoge- 
neous within the sphere of radius R. For a given con- 
figuration of disorder, the LDOS p(rs,w) can be com- 
puted using the dyadic Green function that describes 
the response at point r to the dipole source through 
the relation E(r) = /io w^ G(r,rs,a;)p. It is given by 
p{rs,uj) = 2a;/(7rc^) Im [TrG(rs,rs,aj)], where Tr de- 
notes the trace of a tensor and c is the speed of light in 
vacuum [2l|. In terms of the LDOS, the power radiated 
outside the system, and averaged over the orientation 
of the dipole source, reads P = 7ra;^/(12eo) |ppp(rs,a;). 
Denoting by po and Pq the LDOS and the radiated power 
in vacuum, respectively, we obtain the simple equality 



P(rs,^) ^ P_ 
Po {(^) Po ' 



(1) 



As a characterization of the far-field speckle produced 



by the point source, we consider the angular intensity 
correlation function 



J^>)£(uO^_ 
^("'")-(/(u))(/(uO) ' 



(2) 



where / (u) is the directional radiated power in direction 
u, such that J^^I{u)du = P. From Eqs. ^ and (O, 
the fluctuations of the normalized LDOS can be written 






j^ 11 {I (u) I iW))dudW 



= JS2jJ (I (")> {I ("')> [1 + C (u, u')] dudu'. (3) 

Using the hypothesis of statistical isotropy, the aver- 
aged directional radiated power reduces to (/(u)) — 
(Po/47r) {p (ro, uj)) /pq (w), and the angular intensity cor- 
relation becomes a function of x = u • u'. This allows us 
to simplify Eq. ^ into 



{p{r„uj)f 



= 1 
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C(u-u')dudu' 



C (x) dx. 



(4) 



The correlation function can be expanded on the ba- 
sis of Legendre polynomials in the form C (x) — 
J2n=o'^nPn (x) ■ Siucc Pa{x) = 1, the first term is con- 
stant and corresponds to an infinite range correlation. 
To be consistent with initial considerations of the Co 
correlation [a, [Sj , we define the Cq contribution as that 
given by the constant term, such that oq — Cq. The 
integral in Eq. (|4]) is performed by writing C (x) — 
J2'^=Q o,nPo (x) Pn (x) and using the orthogonality con- 
dition of Legendre polynomials J.^ P„ (x) Pm (x) dx — 
26nm/ (2n. + 1) where Snm is the Kronecker delta. We 
finally obtain 



Co = 



{p{rs,uj)f 



-1 = 



Va.r[p{vs,uj)] 
{p{rs,uj)f 



(5) 



Equation ([5]) shows that the Co speckle correlation and 
the normalized variance of the LDOS at the position of 
the emitter are the same, a result that was first derived in 
Ref.|8| based on a diagrammatic approach. Our derivation 
relies only on considerations of energy conservation, and 
is exact within the only assumption of a non-absorbing 
and statistically isotropic medium. In particular, Eq. ([5]) 
holds in all wave transport regimes, from weakly scatter- 
ing to strongly scattering, including Anderson localiza- 
tion. Another feature of our derivation is that it leads 
to a simple interpretation of the Co correlation. For a 
classical point dipole source, changes in the LDOS cor- 
respond to changes in the power transferred to the en- 
vironment (for a quantum emitter, this corresponds to a 
change in the spontaneous decay rate) . In a non absorb- 
ing medium, this power coincides with the total power 



radiated in the far field, whose angular (or spatial) vari- 
ations generate the speckle pattern. Therefore LDOS 
fluctuations are transferred into global fluctuations of the 
speckle patterns, which themselves are encoded into an- 
gular (or spatial) correlations. 

We have carried out numerical simulations in order 
to analyze the dependance of Co on the local environ- 
ment of the emitter. The scattering medium is a three- 
dimensional cluster of N resonant point scatterers ran- 
domly distributed inside a sphere with radius R. A dipole 
emitter is placed at the center of the cluster (position Ts), 
and is surrounded by a small exclusion volume with ra- 
dius -Ro • The geometry is shown in the inset in Fig. [1] In 
the generation of the random configurations of the disor- 
der (i.e. of the positions of the N scatterers), a minimum 
distance rfo is forced between the scatterers. This permits 
to induce a degree of correlation in the disorder, since this 
amounts to simulating an effective hard-sphere potential 
between scatterers. One can define an effective volume 
fraction f — N {do/2) / {R^ — Rq), that can be taken as 
a measure of the degree of correlation of the disorder (/ 
will be denoted by "correlation parameter" in the follow- 
ing). The scatterers are described by their electric po- 
larizability a (w) = —inc^'y/ [oj^ {cj — cjq + 17/2)] where 
Wo is the resonance frequency, and 7 is the linewidth. 
This corresponds to the polarizability of a resonant non- 
absorbing point scatterer (which is similar to that of 
a two-level atom far from saturation). The scattering 
cross-section is as{uj) = (fc*/67r)|Q;(a;)p, with k = w/c. 
The parameter k£^ measures the scattering strength, 
where £-q = [pas{uj)]~^ is the independent-scattering (or 
Boltzmann) mean-free path, p — N/V being the den- 
sity of scatterers. In all the numerical computations that 
follows, we have taken the following set of parameters: 
N = 100, wo = 3 X IQi^Hz, UJ - LOO = 1 x 10^ Hz, 
7 == 1 X lO^Hz, i? = 1.2 nm and Rq = 0.05 ^im. The 
Boltzmann scattering mean- free path is £b ~ 1-9 |J.m so 
that /c£b = 19 and R/Eb = 0.63. 

The calculation of the statistical distribution of the 
LDOS, from which Co can be deduced, amounts to cal- 
culating the Green function G(r,rs,a;) = Go(r, r^jw) -|- 
S(r, rs,uj) for an ensemble of realizations of the scatter- 
ing medium. Since the free-space Green function Go is 
known analytically, we only need to calculate the Green 
function S(r, rs,wo) which corresponds to the scattered 
field. To proceed, we perform a coupled-dipole numeri- 
cal computation. The field exciting scatterer number j 
is given by the contribution of the dipole source and of 
all other scatterers, leading to a set of 3iV self-consistent 
equations [22| : 

N 

Ej ^ pooJ^Go (rj, r,) p -I- a (w) /c^ ^ Go (r^, r^) E^ (6) 

fe=i 

where Tj is the position of scatterer number j and the 
dependence of the Green functions on u have been omit- 



ted. This linear system is solved numerically for each 
configuration of the disordered medium. Once the excit- 
ing electric field on each scatterer is known, it is possible 
to compute the scattered field at the source position r^ 
and deduce the Green dyadic, from which the LDOS p is 
readily obtained. In this numerical approach, near field 
and far field dipole-dipole interactions and multiple scat- 
tering are taken into account rigorously. 




FIG. 1: (Color online) Statistical distribution of the normal- 
ized LDOS p{rs,u})/ paiuj) for a uncorrelated system with cor- 
relation parameter (effective volume fraction) / = 0.001% in 
a double logarithmic plot. Others parameters are: A'^ = 100, 
cjo = 3 X 10^^ Hz, tj - cjo = 1 X 10® Hz, 7 = 1 X lO'' Hz, 
R — 1.2 [j.m and Ro = 0.05 [an. The calculation are per- 
formed with 3 X 10 configurations. This large number of 
configurations is necessary to correctly describe the tail of 
the distribution. Inset: Schematic view of the system. 

We show in Fig.[T]the statistical distribution of the nor- 
malized LDOS p{rs,ij-')/po{uj) in the case of a correlation 
parameter / = 0.001 % (minimum interparticle distance 
do = 7.5 nm) which corresponds to an almost uncorre- 
lated medium. The curve exhibits a broad distribution, 
with values of p/pa ranging from 0.2 to 1000. This corre- 
sponds to much larger fluctuations than those observed 
in the single scattering regime, where p/po slightly de- 
viates from unity [IJ] . The analysis of the lineshape al- 
lows us to distinguish three regimes. Firstly, the curve 
covers a zone corresponding to p/po < 1, which means 
that some configurations lead to a reduction of the LDOS 
compared to that in free space. This effect has been 
analyzed previously and is due to collective interactions 
in the multiple scattering regime [20[. Recent measure- 
ments of the fluorescence lifetime of emitters at the sur- 
face of a volume scattering disordered medium seem to 
provide evidence of this regime [llj. Secondly, in the 
region p/po > 1, a power- law decay is observed, with 
a statistical distribution behaving as P{p) oc p~^^^ (the 
power law is indicated by the dashed line in Fig. [iJ. As 
originally discussed in Ref. [l5l, this behavior is due to 
near-field interactions between the emitter and the near- 



est scatterer, in the regime of dipole-dipole interactions, 
and in the absence of absorption. Evidence of this regime, 
and of its influence on the Cq correlation, have been re- 
ported recently [U]. Thirdly, in the region p/po ^ 1, 
the tail of the distribution deviates from the power law 
^-3/2 rpj^g onset of the deviation depends on the exclu- 
sion volume with radius Rq that encloses the emitter. In- 
deed, in the single-scattering limit (the emitter interacts 
chiefly with the nearest scatterer), one would observe a 
sharp cut-off indicated by the vertical solid line in Fig. [T] 
This corresponds to the value of p/ po given by the near- 
fleld dipole-dipole interaction between the emitter and 
one scatterer located at a distance Rq {p/po = 62 in the 
present case) [15[ . The observation of a tail beyond this 
single-scattering cut-off is the evidence of near-field in- 
teractions with more than one scatterer. As we shall see 
below, this tail also contains the information on the lo- 
cal environment of the emitter, and in particular on the 
degree of correlation of disorder. 
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FIG. 2: (Color online) Same as Fig. [T] with different values 
of the correlation parameter /. As a result of near-field in- 
teractions, the tail of the distribution is the signature of the 
local environment of the emitter. 

We show in Fig. [5] the statistical distribution of p/po 
for different values of the correlation parameter, ranging 
from 1% to 4.2% (i.e. do ranging from 111 to 180 nm). 
Increasing / amounts to increasing the level of correlation 
in the positions of the scatterers (one imposes an effec- 
tive hard-sphere potential with a longer range) . The tail 
of the distribution is substantially affected by the level 
of correlations in the system. Also note that the part of 
the distribution corresponding to p/po smaller than the 
single-scatterer cut-off remains unchanged. This means 
that the sensitivity of Co to the local environment of the 
emitter is driven by the near-field interactions with the 
surrounding scatterers, this information being encoded 
in the tail of the statistical distribution of the LDOS. Al- 
though this tail corresponds to events with a low proba- 



bility, it is at the core of the Co correlation concept. 

In order to visuahze the influence of the correlation of 
disorder directly on Cq , we have plotted in Fig. [3] the 
values of Cq versus the correlation parameter / obtained 
from numerical simulations (for the same system as in 
Figs. [Hand [2|)- A sharp transition is visible at / ~ 2%, 
the value of Co dropping by a factor of 2. In order to give 
a physical interpretation of this behavior, we have also 
plotted in Fig. [3] the values of Co computed by consider- 
ing the interaction with the nearest scatterer only (black 
dashed curve), and with the two nearest scatterers (blue 
line with markers). For large /, the emitter essentially in- 
teracts with one particle (the red solid line and the black 
dashed curve have a similar behavior) and the Cq cor- 
relation can be understood in simple terms. This is the 
regime found experimentally in Ref. 12. We stress here 
that this regime results from a near-field interaction, so 
that the value of Co depends on local microscopic param- 
eters (it cannot be described with the scattering or trans- 
port mean free path as a single parameter). For small /, 
the probability of getting more than one scatterer in the 
vicinity of the emitter becomes non-negligible, and the 
behavior of Co cannot be explained (even qualitatively) 
with a single-scattering model. One sees that by includ- 
ing the interaction with the two nearest scatterers (blue 
curve with markers), one reproduces nicely the behavior 
of the transition. This result demonstrates the high sen- 
sitivity of Co to the level of correlation of the disorder, 
and the key role of near-field interactions. These inter- 
actions are responsible for the non- universality of Co- 
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FIG. 3: (Color online) Co speckle correlation versus the corre- 
lation parameter /. Red solid line: Full numerical simulation 
with the same parameters as in Fig. [l] Blue solid line with 
markers: Calculation considering only the two nearest scat- 
terers (double scattering). Black dashed line: Calculation 
considering only the nearest scatterer (single scattering). 

In summary, we have have derived the relation between 
the Co speckle correlation and the LDOS fluctuations us- 



ing arguments of energy conservation. This simple and 
exact derivation leads for the first time to an interpreta- 
tion of Co based on the fluctuations of the energy deliv- 
ered by a classical dipole source to a disordered environ- 
ment. Using exact numerical simulations, we have shown 
that Co is essentially a correlation resulting from near- 
field interactions. These interactions give Co its non- 
universal character, that is reflected in its high sensitivity 
to the level of correlation of disorder. 
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